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Abstract 

This  article  provides  a  review  of  brain  tissue  alterations  that  may  be  detectable  using  diffusion 
magnetic  resonance  imaging  MRI  (dMRI)  approaches  and  an  overview  and  perspective  on  the 
modern  dMRI  toolkits  for  characterizing  alterations  that  follow  traumatic  brain  injury  (TBI).  Nonin- 
vasive  imaging  is  a  cornerstone  of  clinical  treatment  of  TBI  and  has  become  increasingly  used  for 
preclinical  and  basic  research  studies.  In  particular,  quantitative  MRI  methods  have  the  potential  to 
distinguish  and  evaluate  the  complex  collection  of  neurobiological  responses  to  TBI  arising  from 
pathology,  neuroprotection,  and  recovery.  dMRI  provides  unique  information  about  the  physical 
environment  in  tissue  and  can  be  used  to  probe  physiological,  architectural,  and  microstructural 
features.  Although  well-established  approaches  such  as  diffusion  tensor  imaging  are  known  to  be 
highly  sensitive  to  changes  in  the  tissue  environment,  more  advanced  dMRI  techniques  have  been 
developed  that  may  offer  increased  specificity  or  new  information  for  describing  abnormalities. 
These  tools  are  promising,  but  incompletely  understood  in  the  context  of  TBI.  Furthermore,  model 
dependencies  and  relative  limitations  may  impact  the  implementation  of  these  approaches  and  the 
interpretation  of  abnormalities  in  their  metrics.  The  objective  of  this  paper  is  to  present  a  basic 
review  and  comparison  across  dMRI  methods  as  they  pertain  to  the  detection  of  the  most  com¬ 
monly  observed  tissue  and  cellular  alterations  following  TBI. 
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1  I  INTRODUCTION 

Despite  the  long  history  of  traumatic  brain  injury  (TBI)  as  a  prevalent 
cause  of  death  and  disability  in  humans,  defining  the  neurobiological 
underpinnings  of  damage  and  recovery  following  TBI  remains  a  central 
challenge.  The  complex  collection  of  physiological,  cellular,  and  molecu¬ 
lar  changes  that  follow  TBI  can  appear  to  be  remarkably  heterogene¬ 
ous,  but  at  the  same  time  they  are  highly  organized  into  coordinated 
responses  such  as  neurodegeneration,  inflammation,  and  regeneration. 
The  corpus  of  histological  studies  spanning  a  variety  of  experimental 
animal  models  of  TBI  have  provided  crucial  insights  about  the  patho- 
mechanisms  and  cellular  alterations  that  accompany  posttraumatic 


tissue  change,  but  considerable  work  remains  to  determine  the  spatio- 
temporal  evolution  of  abnormalities,  interrelationships  among  different 
tissue  responses,  and  their  impact  on  health  and  behavioral  outcomes. 
Noninvasive  imaging  in  animal  models  has  the  potential  to  build  on 
what  is  known  from  histology  by  providing  longitudinal  and  whole- 
brain  information,  but  for  this  approach  to  be  successful  it  is  essential 
to  first  improve  the  understanding  of  how  imaging  abnormalities  corre¬ 
spond  to  tissue  and  cellular  changes. 

Diffusion  magnetic  resonance  imaging  (dMRI)  methods  are  partic¬ 
ularly  promising  for  the  development  of  imaging  markers  of  TBI 
pathology  because  they  are  sensitive  to  microscale  water  displace¬ 
ment  as  a  proxy  for  tissue  environment  geometry  and  provide  a  range 
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Significance 

Noninvasive  magnetic  resonance  imaging  may  enable  important 
observations  about  how  the  brain  changes  following  injury. 
Diffusion  magnetic  resonance  imaging  (dMRI)  is  particularly 
promising  for  studying  traumatic  brain  injury  (TBI)  because  it 
quantitatively  characterizes  the  physical  tissue  environment, 
which  may  change  greatly  after  injury.  While  a  number  of 
conventional  and  advanced  dMRI  approaches  have  been 
developed  that  may  benefit  the  detection  of  post-TBI  abnormal¬ 
ities,  much  remains  to  be  known  about  the  cellular  underpin¬ 
nings  and  relative  value  of  these  methods.  This  article  reviews 
TBI-related  cellular  alterations  that  influence  the  tissue  environ¬ 
ment  and  compares  common  dMRI  tools  for  detecting  such 
changes. 


of  quantitative  scalar  metrics  across  the  whole  brain.  Furthermore, 
dMRI  may  be  combined  with  other  conventional  or  advanced  mag¬ 
netic  resonance  imaging  (MRI)  methods  such  as  arterial  spin  labeling, 
susceptibility-weighted  imaging,  or  a  variety  of  contrast  agent  MRI 
approaches  to  provide  complementary  and  comprehensive  outcome 
measures.  Standard  dMRI  methods  and  especially  diffusion  tensor 
imaging  (DTI)  have  already  demonstrated  sensitive  detection  of  abnor¬ 
malities  in  a  number  of  experimental  models  of  TBI.  In  the  past  dec¬ 
ade,  multiple  advanced  dMRI  approaches  have  extended  beyond  the 
conventional  models  with  the  goals  of  improving  the  physical  descrip¬ 
tion  of  water  diffusion  (e.g.,  by  modeling  “non-Gaussian”  diffusion)  or 
parameterizing  dMRI  with  respect  to  the  expected  biological  environ¬ 
ment  (e.g.,  by  modeling  cellular  compartments  and/or  fiber  geometry). 
These  new  tools  will  be  valuable  if  they  are  able  to  improve  the  sensi¬ 
tivity  or  specificity  of  dMRI  following  TBI;  however,  we  lack  a  system¬ 
atic  understanding  of  how  dMRI  methods  differ  from  one  another  for 
detecting  and  describing  tissue  alterations. 

A  number  of  excellent  reviews  exist  to  describe  the  current  under¬ 
standing  of  cellular  mechanisms  of  TBI  in  general  (Bramlett  &  Dietrich, 
2015;  Pekna  &  Pekny,  2012)  and  within  particular  areas  of  neurobiology 
including  neurodegeneration  (Johnson,  Stewart,  &  Smith,  2013;  Stoica  & 
Faden,  2010),  inflammation  (Burda,  Bernstein,  &  Sofroniew,  2016;  Ziebell 
&  Morganti-Kossmann,  2010),  and  myelin  changes  (Armstrong,  Mierzwa, 
Marion,  &  Sullivan,  2016),  among  others.  As  well,  several  existing  reviews 
have  been  published  regarding  MRI  and  DTI  to  study  human  TBI  (Brody, 
Mac  Donald,  &  Shimony,  2015;  Duhaime  et  al.,  2010;  Hulkower,  Poliak, 
Rosenbaum,  Zimmerman,  &  Lipton,  2013),  and  recently  a  pertinent  over¬ 
view  and  summary  of  advanced  dMRI  tools  and  their  relevance  to  clinical 
outcomes  was  published  (Douglas  et  al.,  2015).  The  focus  of  the  present 
review  is  to  combine  what  is  known  from  work  in  experimental  models 
of  TBI  about  tissue  and  cellular  alterations  that  may  affect  the  physical 
tissue  environment  with  a  comparative  description  of  the  major  methods 
for  dMRI  that  may  be  differentially  sensitive  to  TBI-related  tissue  change 
alongside  several  important  caveats  for  their  use  and  interpretation.  The 
first  section  provides  a  categorical  summary  of  cellular  response  to 
trauma,  emphasizing  alterations  with  microstructural,  architectural,  or 
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neuroanatomical  manifestations  that  may  give  rise  to  detectable  dMRI 
abnormalities,  including  a  review  of  the  existing  dMRI  studies  in  experi¬ 
mental  TBI  models.  The  second  section  contains  a  comparative  overview 
of  presently  available  dMRI  methods  from  standard  approaches  to 
advanced  techniques.  The  objective  of  this  article  is  to  provide  a  refer¬ 
ence  for  the  current  understanding  of  these  topics  as  well  as  a  perspec¬ 
tive  to  help  guide  selection  of  dMRI  tools  based  on  particular  aspects  of 
TBI  questions. 

2  I  THE  PHYSICAL  TISSUE  ENVIRONMENT 
CHANGES  FOLLOWING  TRAUMA 

dMRI  is  able  to  detect  abnormalities  only  when  the  physical  tissue  envi¬ 
ronment  changes  in  a  way  that  substantially  alters  the  movement  of  water 
within  a  voxel,  such  as  the  loss  or  gain  of  cellular  boundaries  or  by  organi¬ 
zation  or  disorganization  of  oriented  structures.  Because  dMRI  measure¬ 
ments  record  signals  from  the  motion  of  water  molecules  in  small  imaging 
volumes  (voxels),  their  sensitivity  is  related  to  the  magnitude  of  structural 
changes  incurred  by  injury.  This  implies  that  dMRI  is  sensitive  to  abnormal 
cerebral  physiology  (e.g.,  edema,  vascular  injury,  disruption  of  water  home¬ 
ostasis),  changes  in  tissue  composition  (e.g.,  increased  or  decreased  cellu- 
larity),  and  alterations  in  cellular  morphology  (e.g.,  glial  reactivity  or  neurite 
density  changes).  On  the  other  hand,  functional  and  molecular  changes 
are  less  likely  to  be  detected  by  dMRI  directly  but  may  be  associated  with 
detectable  changes  indirectly.  Table  1  lists  the  major  categories  of  cellular 
changes  arranged  according  to  cell  type  that  may  give  rise  to  altered  tissue 
environment  and  thus  affect  dMRI  measurements. 

2.1  The  physiologic  response  to  TBI  and  alterations 
of  diffusivity 

Brain  contusion  is  a  common  outcome  of  TBI  and  can  be  readily  identi¬ 
fied  on  T2-weighted  MRI  scans  once  blood  products  have  accumulated 
in  the  parenchyma  or  extracellular  fluid  has  increased  (vasogenic 
edema).  In  the  hyperacute  phase,  dMRI  has  been  shown  to  be  uniquely 
sensitive  to  early  pathophysiologies  that  are  not  detectable  using  con¬ 
ventional  MRI  (Smith  et  al.,  1995).  The  early  sensitivity  of  dMRI  to  TBI 
changes  is  similar  to  observations  in  ischemic  stroke,  for  which  diffusiv¬ 
ity  is  known  to  decrease  robustly  during  the  first  few  hours  in  the 
absence  of  other  MRI  changes  (Moseley  et  al.,  1990).  The  first  dMRI 
studies  of  TBI  were  performed  to  characterize  diffusivity  changes  dur¬ 
ing  this  period  in  experimental  models  of  focal  contusion  and  provided 
evidence  for  both  decreased  (Alsop,  Murai,  Detre,  McIntosh,  &  Smith, 
1996;  Stroop  et  al.,  1998;  Unterberg  et  al..  1997)  and  increased  (Han- 
stock,  Faden,  Bendall,  &  Vink,  1994)  diffusivity  in  the  hours  after  injury. 
While  these  findings  may  seem  contradictory,  experimental  differences 
such  as  injury  severity  can  explain  opposite  changes  in  diffusivity 
(Smith  et  al.,  1995)  and  also  highlight  the  potential  for  dMRI  to  distin¬ 
guish  different  pathophysiological  features  of  tissue.  In  particular, 
mechanisms  of  cellular  damage  including  metabolic  disruption,  beading, 
and  cytotoxic  edema  have  been  proposed  to  explain  acutely  reduced 
diffusivity  following  stroke  and  injury,  while  vasogenic  edema  is  gener¬ 
ally  accepted  to  underlie  increased  diffusivity.  While  cellular  disruption 
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TABLE  1  Cellular  Alterations  Known  to  Follow  TBI  and  Associated  dMRI  Changes 


Cell  type  or 
compartment 

TBI-related  alterations 

Tissue  environment 

Expected  diffusion 
changes 

Major  citations 

dMRI  evidence 

Neurons 

cell  loss  necrosis  and 

atrophy,  cavitation, 

decreased  diffusivity 

Sato  et  al.,  2001; 

Assaf  et  al.,  1997; 

apoptosis 

unmasking 

and  anisotropy, 

increased 

anisotropy 

Chen  et  al.,  2003; 

Coleman,  2005; 

Stoica  &  Faden,  2010 

Van  Putten  et  al.,  2005; 
Immonen  et  al.,  2009; 
Laitinen  et  al.,  2015 

axonal  injury 

axon  morphology 
changes  including 
beading  and 
varicosities 

reduction  in  anisotropy 
and  reduction  in 
diffusion,  especially  in 
the  axial  direction 

Johnson  et  al.,  2013 

Mac  Donald,  Dikranian, 
Bayly,  et  al.,  2007; 

Budde  et  al.,  2009; 

Jiang  et  al.,  2011; 

Li  et  al.,  2011; 

Bennett  et  al.,  2012; 
van  de  Looij  et  al.,  2012 

neural  plasticity 

increased  number  of 

increased  anisotropy 

Bach-y-Rita,  2003; 

Kharatishvili  et  al.,  2007; 

sprouting,  arborization 

coherent  processes  and 
new  collaterals 

and/or  changed 
orientation 

Yiu  &  He,  2006; 

Werner  &  Stevens,  2015; 
Meaney  &  Smith,  2015 

Hutchinson  et  al.,  2012; 
Sierra  et  al.,  2015 

Oligodendrocytes 

demyelination  direct 
damage,  chronic 
pathology 

myelination  repair 
remyelination 

degenerating  or  lost 

regenerating 

decreased  anisotropy 

normalized  anisotropy 

Armstrong  et  al.,  2016 

Jiang  et  at,  2011; 

Budde  et  al.,  2011; 

Li  et  al.,  2014; 

Mac  Donald,  Dikranian, 
Song,  et  al.,  2007 

Astrocytes 

Hypertrophy 

proliferation 
glial  scaring 

increased  number  or 
thickness  of  glial 
processes,  possibly 
organized  or  directional 

increased  cellularity 
dense  glia,  increased 
organization 

increased  or  decreased 
anisotropy,  decreased 
diffusivity 

decreased  diffusivity 
decreased  diffusivity, 
increased  anisotropy 

Norton  et  al.,  1992; 
Sofroniew,  2005,  2009; 
Sofroniew  &  Vinters, 

2010;  Wilhelmsson 
et  al.,  2006;  Suzuki 
et  al.,  2012;  Sun  & 

Jakobs,  2012;  Pekny 
et  al.,  2014;  Burda  et  al., 
2016 

Budde  et  al.,  2011; 

Zhuo  et  al.,  2012; 

Mac  Donald, 

Dikranian,  Song,  et  al., 
2007 

Microglia 

phagocytosis 

neural  repair  and 
support 

amoeboid  stage 
microglia 

rod-microglia 

increased  diffusivity 

possible  increased 
anisotropy 

Kreutzberg,  1996; 

Graeber,  2010; 

Wake  &  Fields,  2011; 
Ziebell  et  al.,  2012; 

Roth  et  al.,  2014 

Intracellular  space 

cytotoxic  edema 

cell  swelling 

decreased  diffusivity 

Moseley  et  al.,  1990; 
Pierpaoli  et  al.,  1993; 
Bramlett  &  Dietrich,  2004 

Hanstock  et  al.,  1994; 
Smith  et  al.,  1995; 

Alsop  et  al.,  1996; 
Unterberg  et  al.,  1997; 
Assaf  et  al.,  1997; 

Stroop  et  al.,  1998; 
Albensi  et  al.,  2000; 

Van  Putten  et  al.,  2005; 
Immonen  et  al.,  2009; 
Frey  et  al.,  2014 

Extracellular  space  vasogenic  edema 

excess  extracellular  fluid 

increased  diffusivity 

Note:  For  each  major  cell  type  of  the  brain  (column  1),  different  types  of  abnormalities  that  have  been  observed  to  follow  experimental  brain  injury  are 
categorized  (column  2).  The  resulting  changes  to  the  tissue  environment  (column  3)  and  water  diffusivity  and  anisotropy  (column  4)  are  given  as  well  as 
relevant  citations  for  the  neurobiological  phenomenon  (column  5)  and  dMRI  observations  (column  6). 


generally  precedes  vasogenic  edema  in  stroke,  this  may  not  always  be 
the  case  in  brain  injury  (Bramlett  &  Dietrich,  2004).  For  instance,  obser¬ 
vations  in  a  model  of  focal  cortical  ischemia  (Pierpaoli  et  al.,  1993)  have 
demonstrated  increased  diffusivity  in  nonischemic  brain  regions  with 
edema  adjacent  to  regions  of  decreased  diffusivity  where  ischemic 


damage  was  later  confirmed  by  histology.  The  implication  of  this  for 
TBI  research  is  that  acutely  increased  diffusivity  may  indicate  brain 
regions  that  undergo  edema  without  cellular  disruption,  and  possibly 
these  areas  will  not  progress  to  degenerative  outcomes,  while  regions 
with  acutely  decreased  diffusivity  are  more  likely  to  have  metabolic  or 
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other  cellular  disruption  that  will  result  in  degeneration.  It  is  important 
to  also  note  that  the  detection  of  these  changes  is  unlikely  to  be  the 
same  for  ex  vivo  dMRI  measurements  as  for  in  vivo  measurements 
given  the  absence  of  physiologic  mechanisms  and  water  homeostasis  as 
well  as  the  reduction  in  extracellular  space  upon  tissue  fixation.  Figure 
1  demonstrates  several  of  these  points  by  showing  DTI  maps  for  in- 
vivo  and  ex-vivo  mouse  and  ferret  brains  following  experimental  TBI. 

Diffusivity  findings  from  these  early  studies  have  been  replicated 
many  times  by  more  advanced  work  and  extended  to  describe  inc¬ 
reased  diffusivity  at  chronic  time  points  corresponding  to  histological 
evidence  of  lesioned  tissue  in  regions  with  abnormal  initial  diffusivity 
(Albensi  et  al.,  2000;  Assaf,  Beit-Yannai,  Shohami,  Berman,  &  Cohen, 
1997;  Van  Putten,  Bouwhuis,  Muizelaar,  Lyeth,  &  Berman,  2005).  In 
addition,  abnormalities  of  diffusivity  following  injury  have  been  associ¬ 
ated  with  meaningful  outcomes  including  atrophy,  seizure  susceptibility 
(Frey  et  al.,  2014),  and  learning  impairment  (Immonen  et  al.,  2009). 

2.2  I  Cellular  alterations  after  TBI  can  affect  dMRI 
scalar  parameters 

We  have  discussed  above  the  ability  of  physiologic  changes  to  influ¬ 
ence  the  magnitude  of  water  diffusion  in  ways  that  are  detectable  by 
basic  dMRI  methods.  With  the  advent  of  DTI  and  other  subsequent 
advanced  dMRI  approaches,  it  has  become  possible  to  probe  changes 
in  the  tissue  environment  arising  from  the  altered  microstructure  and 
local  architecture  of  neurons  and  glia.  This  capability  has  compelled 
several  dozen  studies  of  experimental  TBI  intended  to  identify  and 
characterize  imaging  markers  that  are  associated  with  one  or  more  of 
the  robust  cellular  alterations  that  follow  brain  injury.  Below  are  dis¬ 
cussed  the  most  evident  neuronal  and  glial  alterations  that  may  con¬ 
tribute  to  changes  detectable  using  DTI  or  other  dMRI  approaches. 

2.2.1  |  Neuronal,  axonal,  and  myelin  alterations 

Neuron  cell  bodies,  dendrites,  and  axons  encompass  a  wide  range  of 
spatial  dimensions  and  morphologic  features  that  may  affect  the  char¬ 
acteristics  of  water  diffusion  in  tissue.  For  example,  hippocampal  sub¬ 
fields  can  be  distinguished  on  the  basis  of  different  dMRI  values 
observed  to  be  related  to  dendrite  morphology  (Shepherd,  Ozarslan, 
King,  Mareci,  &  Blackband,  2006),  and  cortical  anisotropy  during  devel¬ 
opment  has  been  shown  to  reflect  changes  in  neurite  density  and  com¬ 
plexity  (Huang,  Yamamoto,  Hossain,  Younes,  &  Mori,  2008).  While 
neuronal  features  in  gray  matter  tissue  regions  may  be  detectable  using 
dMRI,  the  far  more  common  use  is  to  characterize  the  myelinated 
axons  of  the  brain’s  white  matter  pathways  using  diffusion  anisotropy 
measures  that  are  highly  sensitive  to  cylindrical  geometry.  Following 
TBI,  neuron  loss,  axonal  damage,  demyelination,  and  neural  reorganiza¬ 
tion  may  drastically  alter  tissue  composition  as  well  as  the  local  tissue 
geometry  in  ways  that  are  detectable  using  dMRI. 

Neuronal  cell  death 

The  irreversible  loss  of  neurons  defines  neurodegenerative  disorders 
resulting  in  the  loss  of  brain  function,  and  TBI  has  been  associated  with 
several  mechanisms  of  neuronal  death,  each  affecting  the  time  course 
and  tissue  environment  of  the  brain  following  injury  (Chen,  Pickard,  & 
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Harris,  2003;  Sato,  Chang,  Igarashi,  &  Noble,  2001;  Stoica  &  Faden, 
2010).  Neuronal  death  may  result  from  direct  immediate  physical  dam¬ 
age  induced  by  biomechanical  forces  that  disrupt  cell  membranes,  caus¬ 
ing  primary  cell  loss,  or  may  be  initiated  by  a  range  of  secondary  cell  loss 
processes  that  lead  to  active  or  programmed  cell  death.  The  loss  of 
neurons  in  the  tissue  environment  should  ultimately  lead  to  increased 
diffusivity  from  the  loss  of  cell  membranes  and  gross  morphological 
alterations  such  as  cavitation  and  atrophy,  but  cell  death  mechanisms 
that  precede  these  changes  may  or  may  not  encompass  morphological 
features,  such  as  cell  swelling,  that  affect  diffusivity  (Stoica  &  Faden, 
2010).  When  neuron  death  occurs,  the  axon  and  myelin  sheath  degen¬ 
erate  rapidly  in  a  process  known  as  Wallerian  degeneration,  which  can 
lead  to  atrophy  in  the  chronic  phase.  In  the  hours  after  injury,  an  accu¬ 
mulation  of  organelles  in  the  axon  causes  swelling,  and  the  myelin 
sheath  is  stretched.  This  is  followed  by  disruption  of  the  cytoskeleton 
and  granular  disintegration  that  completely  degrades  the  axon  within 
several  days  following  injury.  Morphological  changes  that  arise  from 
Wallerian  degeneration  and  that  may  influence  the  tissue  environment 
are  in  common  with  those  described  below  for  direct  axonal  damage 
and  myelin  pathology.  In  fact,  there  is  evidence  and  strong  support  for 
a  convergent  set  of  degeneration  mechanisms  following  a  diverse 
range  of  initial  neuron  injury  (Coleman,  2005). 

Evidence  for  the  contribution  of  cell  loss  to  dMRI  changes  has 
been  demonstrated  in  a  number  of  studies  by  comparing  dMRI  metrics 
with  cellularity  (Van  Putten  et  al.,  2005)  or  by  correlation  of  the  MRI 
and  histologically  defined  values  for  lesion  volume  (Assaf  et  al.,  1997). 
For  robust  dMRI  changes  in  the  acute  stage,  numerous  tissue  altera¬ 
tions  are  likely  present  in  these  regions  of  primary  damage  and  cellular 
loss  near  the  injury  site.  The  impact  of  secondary  cell  loss  on  the  tissue 
environment  during  the  chronic  phase  has  been  documented  and  asso¬ 
ciated  with  dMRI  outcomes  for  several  regions  including  the  increased 
mean  diffusivity  (MD)  in  the  hippocampus  (Laitinen,  Sierra,  Bolkvadze, 
Pitkanen,  &  Grohn,  2015)  and  increased  fractional  anisotropy  (FA)  in 
the  thalamus  (Immonen  et  al.,  2009). 

Traumatic  and  diffuse  axonal  injury 

Biomechanical  forces  in  the  brain  caused  by  trauma  preferentially 
injure  the  long  thin  axonal  processes  of  neurons  that  connect  different 
neuroanatomical  regions  (Johnson  et  al.,  2013).  While  some  axonal  loss 
may  result  from  direct  disruption  of  the  cell  membrane,  the  majority  of 
axonal  damage  is  secondary  and  follows  a  course  of  cellular  changes 
that  are  most  robust  in  the  days  after  injury,  but  have  been  found  to 
persist  for  much  longer  and  may  lead  to  axonal  degeneration  or  func¬ 
tional  impairment.  While  gross  degeneration  of  axons  may  cause 
changes  in  the  composition  of  tissue,  typically  a  subset  of  axons  within 
a  tract  are  affected.  The  effects  of  axonal  damage  that  are  most  conse¬ 
quential  for  detection  by  dMRI  are  changes  in  the  axonal  geometry  to 
which  diffusion  techniques  are  remarkably  sensitive.  In  particular,  dis¬ 
ruption  of  axonal  transport  from  cytoskeletal  stretching  leads  to  the 
appearance  of  axonal  varicosities  and  retraction  bulbs  along  the  axon, 
which  can  reduce  diffusivity  especially  along  the  axon  and  increase  the 
amount  of  restricted  water  contributing  to  the  diffusion  signal.  Early 
DTI  studies  of  mouse  white  matter  identified  decreases  in  anisotropy 
and  diffusivity  along  the  axon  (axial  diffusivity,  Dax)  4  to  24  hr 
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FIGURE  1  Two  examples  of  DTI  metric  abnormalities  following  experimental  TBI  in  ferret  (a)  and  mouse  (b)  brains.  For  each  species,  the 
in  vivo  and  ex  vivo  FA  and  trace  (TR)  maps  and  T2-weighted  images  are  shown  from  the  same  animal  after  controlled  cortical  impact  (CCI 
site  indicated  by  red  arrowhead).  Several  key  features  of  diffusion  changes  after  TBI  are  demonstrated  in  this  figure  including  heterogeneity 
of  diffusivity  abnormalities  within  regions  of  edema  shown  by  values  of  TR  that  are  increased  (a),  decreased  (b),  or  normal  (a  and  b)  within 
tissue  regions  with  T2  hyperintensity.  Distinct  profiles  of  TR  and  FA  can  also  be  found  in  this  figure  by  comparing  images  in  the  middle  row 
where  TR  is  relatively  normal  for  both  the  ferret  and  mouse  brains,  but  FA  is  decreased  in  the  ferret  brain  white  matter  at  1  week  (a)  and 
increased  in  the  mouse  brain  cortex  at  12  weeks  (b).  By  comparing  the  middle  and  last  rows  of  in  vivo  and  ex  vivo  maps  from  the  same  ani¬ 
mal  at  the  same  time  point,  a  distinct  pattern  can  be  found  for  the  ferret  brain  at  1  week  (a)  in  which  subdomains  of  increased  TR  (near  the 
red  arrowhead)  and  decreased  TR  (yellow  arrow)  can  be  found  in  regions  of  unremarkable  in  vivo  TR.  In  contrast,  the  same  region  of 
increased  FA  can  be  found  in  both  the  in  vivo  and  ex  vivo  mouse  brain  12  weeks  after  CCI  (b).  The  observations  depicted  in  this  figure 
demonstrate  several  of  the  key  points  described  in  the  text 


following  controlled  cortical  impact  (CCI)  consistent  with  traumatic 
axonal  injury  pathology  and  confirmed  by  amyloid  beta  precursor  pro¬ 
tein  and  neurofilament  staining  (Mac  Donald,  Dikranian,  Song,  et  al., 
2007)  in  the  absence  of  major  changes  in  myelin  basic  protein  expres¬ 
sion.  This  is  consistent  with  findings  in  a  study  of  isolated  axonal  injury 
found  to  be  related  to  decreased  FA  and  axial  diffusivity  (Budde,  Xie, 
Cross,  &  Song,  2009).  In  addition,  numerous  subsequent  experimental 
TBI  studies  of  CCI  (Budde,  Janes,  Gold,  Turtzo,  &  Frank,  2011; 
Davoodi-Bojd  et  al.,  2014;  Jiang  et  al.,  2011;  Li  et  al.,  2014;  Wang 
et  al.,  2013;  Zhuo  et  al.,  2012),  closed  head  injury  (CHI)  (Bennett,  Mac 
Donald,  &  Brody,  2012;  Li,  Li,  Feng,  &  Gu,  2011;  van  de  Looij  et  al., 
2012),  and  blast  (Budde  et  al.,  2013;  Calabrese  et  al.,  2014)  have  been 
able  to  detect  histologically  verified  axonal  pathology  using  dMRI 
approaches. 

Myelin  damage  and  loss 

Neuronal  degeneration  and  axonal  damage  are  intricately  related  to 
pathophysiology  and  structural  features  of  the  myelin  sheath 
(Armstrong  et  al.,  2016),  which  has  considerable  influence  on  the  dMRI 
signal  (Beaulieu,  2002).  If  a  neuron  is  damaged  to  the  point  of  discon¬ 
nection  and  degeneration,  the  myelin  sheath  collapses  irreversibly  and 
the  remaining  debris  lingers  in  the  interstitial  space  until  it  is  eventually 
cleared.  Demyelination  is  a  different  process  by  which  direct  damage 
causes  apoptosis  of  the  oligodendrocytes  that  form  the  myelin  sheath, 
and  often  a  subset  of  demyelinated  axons  are  interspersed  within  a 
white  matter  tract.  Unlike  axonal  degeneration,  this  process  can  be 
reversible  with  remyelination  and  restoration  of  neuronal  function. 

As  described  in  the  previous  section,  TBI-related  changes  to 
myelinated  axons  are  unlikely  to  present  a  single  pathologic  feature, 
but  the  relative  contribution  of  axonal  and  myelin  alterations  may 


accessible  to  dMRI  methods.  For  example,  while  reduced  axial  diffusiv¬ 
ity  was  reported  to  acutely  follow  CCI  in  the  presence  of  axonal  dam¬ 
age  without  myelin  changes,  a  second  study  in  the  same  model 
(Mac  Donald,  Dikranian,  Bayly,  Holtzman,  &  Brody,  2007)  found  that 
weeks  after  injury,  when  demyelination  became  evident,  reduced  ani¬ 
sotropy  was  no  longer  accompanied  by  reduced  axial  diffusivity  but, 
instead,  by  increased  radial  diffusivity.  This  is  supported  by  similar  find¬ 
ings  of  increased  radial  diffusivity  in  a  model  with  demyelination  with¬ 
out  axonal  damage  (Sun  et  al.,  2006)  as  well  as  subsequent  studies 
demonstrating  DTI  changes  primarily  in  anisotropy  with  abnormal  mye¬ 
lin  staining  (Budde  et  al.,  2011;  Jiang  et  al.,  2011;  Li  et  al.,  2014).  Nota¬ 
bly,  changes  in  myelination  and  myelinated  fibers  may  influence 
anisotropy  differently  depending  on  the  tissue  investigated.  For  exam¬ 
ple,  loss  of  myelinated  fibers  in  white  matter  was  observed  to  reduce 
FA  during  the  chronic  period,  while  their  loss  in  the  thalamus  resulted 
in  increased  FA  (Laitinen  et  al.,  2015). 

Sprouting,  remodeling,  and  regeneration 

The  neuronal  response  to  brain  injury  also  includes  remodeling  of  neuro¬ 
circuitry  and  neurite  plasticity  with  the  potential  to  promote  recovery 
(Bach-y-Rita,  2003;  Werner  &  Stevens,  2015).  These  processes  appear 
to  be  more  complex  and  variable  than  those  related  to  neural  damage 
and  highly  modulated  by  molecular  signaling  in  the  cell  environment  for 
both  promotion  and  inhibition  of  axonal  growth  (Werner  &  Stevens, 
2015;  Yiu  &  He,  2006).  Nevertheless,  there  is  a  clear  behavioral  time 
course  for  functional  recovery  that  implies  the  presence  of  neural  repair 
and/or  circuit  remodeling  in  the  brain.  At  the  cellular  level,  a  time  course 
has  been  proposed  in  which  synaptic  changes  drive  remodeling  in  the 
days  and  weeks  following  injury,  and  more  widespread  connections  in 
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FIGURE  2  Cross-model  comparison  of  scalar  maps  in  the  injured  brain.  A  range  of  tissue  and  injury-related  contrasts  may  be  visually 
observed  in  this  collage  of  16  representative  metrics  in  the  same  slice  from  different  dMRI  models.  This  cross-model  view  of  scalar  maps 
demonstrates  the  potential  for  nonredundant  information  about  regions  of  injury  that  may  be  gleaned  from  different  models.  DTI  metrics  of 
fractional  anisotropy  (FA),  trace  (TR),  axial  and  radial  diffusivity  (Dax  and  Drad),  directionally  encoded  color  (DEC)  map  weighted  by  lattice 
index,  DEC  weighted  by  Westin  linear  anisotropy  (WL)  and  DEC  weighted  by  Westin  planar  anisotropy  (WP),  DKI  metrics  of  mean  kurtosis 
(MK),  axial  and  radial  kurtosis  (AK  and  RK)  and  kurtosis  FA  (KFA),  MAP-MRI  metrics  of  return  to  the  origin,  axis,  and  plane  probabilities 
(RTOP,  RTAP,  and  RTPP),  propagator  anisotropy  (PA)  and  non-Gaussianity  (NG)  and  NODDI  metrics  of  compartment  volume  fractions  for 
isotropic  free  water  (Viso),  intracellular  water  (Vic)  and  intracellular  restricted  water  (Vir),  and  orientation  dispersion  index  (ODI).  Insets  of 
each  map  show  tissue  near  the  injury  site  where  dMRI  values  are  expected  to  be  abnormal 


neurocircuitry  become  evident  during  the  chronic  period  following  injury 
(Meaney  &  Smith,  2015). 

The  most  evident  relationship  between  dMRI  metrics  and  neural 
reorganization  following  TBI  has  been  found  in  the  hippocampus  fol¬ 
lowing  focal  experimental  TBI  in  rats.  In  the  hippocampus,  a  pattern  of 
acutely  decreased  and  chronically  increased  diffusivity  has  been 
reported  (Kharatishvili,  Immonen,  Grohn,  &  Pitkanen,  2007)  as  well  as 
changes  in  anisotropy  and  tissue  orientation  on  the  side  of  injury 
(Hutchinson,  Rutecki,  Alexander,  &  Sutula,  2012;  Sierra,  Laitinen, 


Grohn,  &  Pitkanen,  2015).  In  several  of  these  studies,  the  observed 
dMRI  changes  were  associated  with  mossy  fiber  sprouting,  which  is  a 
well-established  rewiring  of  hippocampal  circuitry  associated  with  epi- 
leptogenesis  (Sutula  et  al.,  1998)  and  known  to  result  following  TBI 
(Golarai,  Greenwood,  Feeney,  &  Connor,  2001).  While  it  is  appealing 
that  increases  in  FA  can  be  used  as  a  marker  of  neuronal  plasticity, 
great  care  must  be  taken  in  making  this  interpretation  as  other  sources 
of  anisotropy  are  present  in  the  gray  matter  that  may  change  with 
injury  such  as  organized  gliosis  (Budde  et  al.,  2011). 
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2.2.2  Glial  alterations 

As  more  becomes  known  about  the  function  of  glial  cells  in  the  central 
nervous  system  (CNS),  it  is  clear  that  they  are  vital  to  a  range  of  func¬ 
tional  processes  in  the  healthy  brain  including  in  homeostatic  regula¬ 
tion,  synapse  structure  and  signaling,  vascular  coupling,  and  waste 
removal,  among  others.  Astrocytes  are  at  least  as  numerous  as  neurons 
in  the  brain  (Herculano-Houzel,  2014;  Sofroniew  &  Vinters,  2010)  with 
nonoverlapping,  but  functionally  connected  astrocyte  domains  of  55 
microns  in  mice  to  145  microns  in  humans  (Oberheim  et  al.,  2009;  Sun 
&  Jakobs,  2012),  and  each  astrocyte  exerts  influence  over  most  of  the 
cellular  elements  in  its  domain,  especially  synapses.  White  matter 
astrocytes  (fibrous)  have  smaller  cell  bodies,  fewer  and  less  branched 
processes,  and  elongated  morphology  along  axons  compared  with  gray 
matter  (protoplasmic)  astrocytes,  and  they  also  stain  more  readily  for 
glial  fibrillary  acidic  protein  (GFAP)  (Lundgaard,  Osorio,  Kress,  Sang- 
gaard,  &  Nedergaard,  2014;  Sofroniew  &  Vinters,  2010).  Microglia  are 
smaller  than  astrocytes,  and  unlike  astrocytes,  they  are  not  directly 
connected  to  neural  or  vascular  networks;  rather,  each  microglial  cell  in 
healthy  tissue  surveys  its  own  territory  with  mobile  processes  that  are 
highly  sensitive  to  changes  in  the  microenvironment.  While  astrocytes, 
microglia,  and  other  “supporting  cells”  of  the  CNS  are  far  from  quies¬ 
cent  in  the  healthy  brain,  their  response  to  injury  is  by  comparison 
robust,  diverse,  and  coordinated  across  multiple  cell  types  with  major 
roles  in  the  consequences  of  injury  (Burda  &  Sofroniew,  2014).  Fur¬ 
thermore,  changes  in  glial  cells  are  remarkably  morphological  in  nature, 
demonstrating  cytoarchitectural  alterations  that  seem  likely  to  influ¬ 
ence  the  microstructural  tissue  environment  in  ways  that  are  detecta¬ 
ble  by  dMRI. 

Astrocyte  reactivity:  Hypertrophy,  proliferation, 
and  glial  scar  formation 

Astrocytes  become  “reactive”  in  response  to  a  range  of  triggering 
events  or  environmental  conditions  including  neuroinflammation,  ische¬ 
mia,  and  mechanical  injury  (Burda  et  al.,  2016;  Sun  and  Jakobs,  2012). 
Reactive  astrocytosis  plays  a  major  role  in  neuroprotection,  especially 
during  the  acute  period,  but  can  also  hamper  regenerative  processes 
and  recovery  from  injury  at  later  time  points  (Pekny,  Wilhelmsson,  & 
Pekna,  2014;  Sofroniew,  2005).  Following  TBI,  astrocyte  reactivity  is 
dependent  on  the  nature  and  severity  of  injury  and  marked  by  hyper¬ 
trophy  and  proliferation  during  the  acute  period  that  may  either  sub¬ 
side  or  lead  to  the  formation  of  glial  scaring  during  the  chronic  period. 
Hypertrophy  of  astrocytes  is  a  hallmark  feature  of  reactivity  driven  by 
the  upregulation  of  intermediate  filament  proteins  (e.g.,  vimentin  and 
GFAP)  and  likely  to  be  consequential  for  dMRI  measurements,  but  the 
interpretation  of  staining  for  GFAP,  a  common  marker  of  glial  reactiv¬ 
ity,  has  several  caveats.  In  particular,  GFAP  staining  is  specific  to  the 
cytoskeleton  and  labels  only  15%  of  the  astrocyte  volume,  so  the  stain¬ 
ing  pattern  may  not  represent  the  true  cellular  morphometry  (Sun  and 
Jakobs,  2012).  Nevertheless,  GFAP  immunohistochemistry  (IHC)  is  an 
important  indicator  of  changes  in  the  number,  length,  and  thickness  of 
astrocyte  processes  to  measure  hypertrophic  morphology  and  has 


been  used  successfully  to  describe  quantitative  changes  in  morphology 
that  follow  experimental  TBI  (Norton,  Aquino,  Hozumi,  Chiu,  & 
Bronson,  1992). 

Proliferation  or  migration  of  astrocytes  near  injury  would  also 
potentially  affect  the  environment  by  changing  tissue  composition. 
Whether  astrocytes  actively  migrate  toward  the  site  of  injury  is  not 
presently  clear  (Wilhelmsson  et  al.,  2006),  but  the  proliferation  of 
new  astrocytes  has  been  demonstrated  following  TBI,  although  it  is 
limited  to  severe  injury  and  typically  localized  near  the  injury  site 
and  involved  in  the  formation  of  a  glial  scar  (Sofroniew,  2009; 
Suzuki,  Sakata,  Kato,  Connor,  &  Morita,  2012).  In  milder  injury, 
another  caveat  for  the  interpretation  of  GFAP  staining  has  been 
observed  when  a  greater  number  of  GFAP-positive  cells  are  found 
after  injury  that  are  not  from  new  astrocytes  but,  rather,  from 
upregulation  of  GFAP  in  cells  that  would  not  normally  stain  positive 
(Sofroniew  &  Vinters,  2010). 

Formation  of  a  glial  scar  is  the  chronic  outcome  of  reactive 
astrocytosis  and  provides  a  physical  barrier  to  cordon  off  damaged 
or  toxic  tissue  to  protect  viable  tissue,  and  while  the  conventional 
view  is  that  the  glial  scar  also  inhibits  regenerative  growth  by  cell 
signaling  (Silver  &  Miller,  2004),  this  paradigm  has  recently  been 
called  into  question  by  evidence  that  it  may  instead  promote  regen¬ 
eration  in  some  cases  (Anderson  et  al.,  2016).  The  glial  scar  itself  is 
a  very  dense  region  of  astrocytes  absent  of  neuronal  projections 
and  vasculature.  In  such  an  environment,  the  diffusion  of  water  is 
expected  to  be  decreased  compared  with  healthy  tissue.  Further¬ 
more,  astrocytes  in  the  glial  scar  typically  demonstrate  a  highly  ori¬ 
ented  and  elongated  morphology,  which  has  been  observed  to 
increase  diffusion  measurements  of  anisotropy  in  experimental  TBI 
(Budde  etal.,2011). 

Several  dMRI  studies  have  demonstrated  a  strong  relationship 
between  dMRI  measures  and  astrocytic  changes  following  trauma 
indicated  by  GFAP  staining  (Budde  et  al.,  2011;  Mac  Donald,  Dikra- 
nian,  Bayly,  et  al.,  2007;  Zhuo  et  al.,  2012).  While  early  decreases  in 
diffusivity  have  been  shown  to  be  colocalized  with  increased  GFAP 
staining  in  cortical  tissue  near  the  injury  site  (Zhuo  et  al.,  2012),  a 
different  profile  of  dMRI  change  has  been  observed  during  the 
chronic  period  in  which  anisotropy  is  increased  and  associated  with 
increased  anisotropy  of  GFAP  staining  in  the  same  region  of  gray 
matter  near  the  injury  site  (Budde  et  al.,  2011).  In  the  white  matter, 
astrocyte  reactivity  was  observed  to  accompany  axonal  damage  and 
associated  with  decreases  in  FA  during  the  chronic  period,  but  it  did 
not  contribute  to  decreased  FA  in  the  acute  period  (Mac  Donald, 
Dikranian,  Bayly,  et  al.,  2007).  Taken  together,  these  studies  demon¬ 
strate  the  heterogeneity  of  DTI  outcomes  that  can  arise  from  the 
same  pathomechanism  depending  on  the  time  after  injury  and  tissue 
type  affected.  Along  with  the  notable  lack  of  specificity  of  dMRI 
metrics— that  is,  multiple  cellular  alterations  could  lead  to  the  same 
DTI  abnormality— the  development  of  advanced  dMRI  methods  that 
are  able  to  disentangle  various  features  of  tissue  change  would  be 
consequential  for  improved  identification  and  interpretation  of 
pathology  following  TBI. 
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TABLE  2  Comparative  Overview  of  Diffusion  MRI  Models 


Acquisition 

Modeling  approach 

DWI  sampling 

DWI  weighting 

Metrics 

Seminal  references 

ADC 

S  =  S0e-bD 

Few  DWIs;  one  b=0  and 
1  DWI 

Low 

ADC 

Eccles  et  al.,  1988 

DTI 

S  =  S0e~b:D 

Few  DWIs;  one  shell, 
minimum  six  directions 

Low 

Combinations  of  11,  12, 
13.  e.g.  FA,  TR,  WL,  WP 

Basser  et  al.,  1994 

DKI 

S=S0e-bD+y°2K 

Moderate  number  of 
DWIs;  two  shells 

Low  and  moderate  only 

DTI  metrics  and  mean 
kurtosis,  axial  and  radial 
kurtosis  and  KFA 

Jensen  et  al.,  2005; 

Tabesh  et  al.,  2011; 

Glenn  et  al.,  2015 

MAP 

Asymmetric  simple 
harmonic  oscillator 
reconstruction  and 
estimation 

Moderate  to  many 

DWIs;  multishell 
acquisition 

Low,  moderate,  and  high 

DTI  metrics  and 
non-Gaussianity, 
zero-displacement 
probabilities,  propagator 
anisotropy,  ODFs 

Ozarslan  et  al.,  2013; 

Avram  et  al.,  2016 

DSI 

Model-free 

Many  DWIs;  Cartesian 
grid 

Low,  moderate,  and  high 

ODFs  possible  to 
generate  zero- 
displacement 
probabilities 

Tuch  et  al.,  2003 

Q-ball 

Model-free 

Moderate  no.  of  DWIs; 
single-shell  HARDI  ac¬ 
quisition 

High 

ODFs  possible  to 
generate  zero- 
displacement 
probabilities 

Tuch,  2004 

CHARMED 

Intra/extra-axonal 
compartments  modeled 
by  restricted/hindered 
sheets  and  cylinders 

Multishell  acquisition 

Low,  moderate,  and  high 

Restricted  and  hindered 
component  fractions; 
cone  of  uncertainty 

Assaf  &  Basser,  2005 

Axcaliber 

Similar  to  CHARMED, 
but  with  additional 
modeling  of  axon 
diameter 

Multishell  acquisition 

Low,  moderate,  and  high 
flexible 

CHARMED  metrics  and 
axon  diameter 

Assaf  et  al.,  2008 

NODDI 

Watson  distributed 
cylinders  and  sticks 

Moderate  no.  of  DWIs: 

multishell 

acquisition 

Low,  moderate,  and  high 
(flexible) 

Cellular  fractions, 
orientation  dispersion 
index 

Zhang  et  al.,  2012; 

Tariq  et  al.,  2016 

WMTI 

Intra/extra-axonal 
compartments  modeled 
with  the  Gaussian  part 
of  the  DKI  model 

Moderate  no.  of  DWIs; 
two  shells  (same  as  DKI) 

Low  and  moderate  only 

Axonal  water  fraction, 
intra-axonal  diffusivity, 
extra-axonal  radial/axial 
diffusivity,  extra-axonal 
tortuosity 

Fieremans  et  al.,  2011 

CSD  tractography  Constrained  spherical 
deconvolution 

Single-shell  HARDI 
acquisition 

High 

ODFs  and  tractograms 

Tournier  et  al.,  2004,  2012 

Note:  For  physical  and  biophysical  dMRI  models  (column  1,  light  and  dark  gray  shading,  respectively),  a  summary  of  the  modeling  approach  (column 
2)  and  acquisition  strategy  (columns  3  and  4)  is  given  as  well  as  the  primary  scalar  metrics  (column  5).  Abbreviations:  ADC  -  Apparent  Diffusion 
Coefficient;  DWI  -  diffusion  weighted  image;  DTI  -  Diffusion  Tensor  Imaging;  TR  -  Trace;  WL  -  Westin's  Linear  Anisotropy;  WP  -  Westin's  Planar 
Anisotropy;  DKI  -  Diffusion  Kurtosis  Imaging;  KFA  -  Kurtosis  FA;  MAP  -  Mean  Apparent  Propagator;  ODF  -  Orientation  Distribution  Function;  DSI 
-  Diffusion  Spectrum  Imaging;  CHARMED  -  Composite  Hindered  And  Restricted  Model  of  Diffusion;  NODDI  -  Neurite  Orientation  Dispersion  and 
Density  Imaging;  WMTI  -  White  Matter  Tract  Integrity;  CSD  -  Constrained  Spherical  Deconvolution;  HARDI  -  High  Angular  Resolution  Diffusion 
Imaging. 


Microglial  reactivity:  Altered  form  and  function 
of  the  immune  cells  of  the  CNS 

Microglia  are  early  responders  to  tissue  environment  change  and  the 
primary  cell  type  for  the  brain’s  immune  response  (Graeber,  2010; 
Kreutzberg,  1996;  Wake  &  Fields,  2011).  Like  astrocytes,  microglia 
become  "reactive”  in  response  to  injury  in  a  graded  manner,  but 
with  morphological  features  that  are  quite  distinct.  The  classic  pat¬ 
tern  of  microglial  reactivity  is  defined  by  morphological  stages  (Zie- 
bell,  Adelson,  &  Lifshitz,  2015)  whereby  microglial  processes  retract 
and  the  cell  body  enlarges  such  that  the  resting  state  or  “ramified” 


microglial  morphology  with  extensive  processes  is  quite  distinct 
from  the  reactive  or  “amoeboid"  state  with  morphology  that  is 
nearly  spherical.  In  addition  to  these  hallmark  phenotypes  of  micro¬ 
glial  activation,  other  observations  of  morphological  alteration  are 
known  to  follow  injury  including  “honeycomb”-  or  “jellyfish”-like 
arrangement  of  microglial  processes  (Roth  et  al.,  2014)  acutely  and 
rod  morphology  that  is  most  evident  a  week  after  injury  (Ziebell, 
Taylor,  Cao,  Harrison,  &  Lifshitz,  2012).  Clearly,  the  morphological 
response  of  microglia  is  remarkably  diverse  and  appears  to  reflect 
information  about  the  tissue  state  following  injury,  making  this  cell 
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FIGURE  3  Two  examples  of  caveats  are  shown  for  the  biophysical  representation  of  diffusion  MRI  information  by  tractography  in  the 
ferret  (a)  and  mouse  (b)  brain  using  the  same  approach  and  parameters.  In  the  ferret  brain  near  the  site  of  a  penetrating  injury,  FA  is  low 
and  few  tracts  can  be  found  in  the  body  of  the  white  matter  compared  with  the  contralateral  side.  However,  investigation  of  this  region 
using  IHC  of  the  same  brain  reveals  the  presence  of  myelinated  axons  (indicated  by  MOG  IHC— myelin  oligodendrocyte  glycoprotein)  and 
upregulated  staining  of  astrocytes  (by  GFAP  IHC— glial  fibrillary  acidic  protein).  The  interpretation  of  the  tractography  in  this  case  could 
indicate  a  loss  of  white  matter  fibers,  when  in  fact  the  underlying  pathology  appears  to  be  more  related  to  gliosis.  In  the  mouse  brain  (b),  a 
region  of  increased  FA  and  aberrant  “tracts"  can  be  found  in  the  cortex  near  the  injury  site;  however,  inspection  by  IHC  reveals  a 
disruption  of  MOG  staining  and  upregulation  and  organized  GFAP  staining  in  this  tissue  region  of  the  same  animal.  The  interpretation  of 
tractography  in  this  case  could  suggest  cortical  plasticity,  when  in  fact  the  underlying  alteration  is  more  related  to  glial  changes.  This  is 
similar  to  a  finding  reported  by  Budde  et  al.  (2011).  Taken  together,  this  figure  emphasizes  the  need  for  careful  interpretation  of  dMRI 
findings,  especially  for  biophysical  models  such  as  tractography,  which  directly  report  neurobiological  metrics  based  model  assumptions 


type  a  promising  target  for  examination  by  dMRI,  although  the  small 
size  of  microglial  cells  may  limit  their  influence  on  the  tissue 
environment. 

3  I  THE  MODERN  DMRI  TOOLKIT  AND 
DETECTION  OF  POST-TBI  TISSUE  CHANGE 

It  is  evident  from  the  previous  section  that  a  substantial  and  growing 
body  of  dMRI  studies  has  emerged  that  demonstrates  the  sensitivity  of 
diffusion  weighted  imaging  (DWI)-  and  DTI-derived  diffusivity  and  ani¬ 
sotropy  metrics  for  detecting  TBI-related  brain  abnormalities  across  a 
wide  range  of  changes  in  the  brain  tissue  environment.  However,  a 
fundamental  limitation  of  these  studies  is  the  low  specificity  of  DTI 
metrics  for  identifying  particular  cellular  features  of  interest.  In  part, 
this  can  be  addressed  by  studies  that  incorporate  histological  outcomes 
to  substantiate  the  interpretation  of  DTI  changes.  Another  promising 
avenue  for  improving  the  utility  of  dMRI  in  TBI  research  is  application 
of  more  advanced  modeling  methods  that  may  provide  additional  or 
more  specific  measures  of  tissue  and  cellular  changes  following  TBI, 
and  a  number  of  potentially  advantageous  dMRI  approaches  have  been 
developed  in  recent  years  (see  Figure  2  for  representative  scalar  maps 
across  different  diffusion  models).  To  this  end,  this  section  presents  a 
concise  and  systematic  overview  of  conventional  and  advanced  dMRI 
methods  and  a  perspective  on  their  relative  strengths  and  limitations  in 
the  context  of  TBI  research. 


The  general  pipeline  that  is  used  to  perform  most  dMRI  techniques 
is  (1)  acquisition  of  DWI  volumes  having  particular  diffusion  weighting 
direction  sets  and  b-values,  (2)  correction  of  DWI  volumes  for  artifacts 
and  distortions,  (3)  modeling  of  the  DWI  data  (or  application  of  a  math¬ 
ematical  transform)  to  generate  scalar  metric  maps,  and  (4)  analysis  and 
interpretation  of  the  scalar  metric  maps  to  answer  an  experimental 
question.  While  this  pipeline  appears  to  be  sequential  in  nature,  there 
is  interdependence  of  each  step  on  the  others  that  requires  careful 
planning  of  all  steps  prior  to  data  collection.  Table  2  outlines  basic 
information  related  to  these  steps  across  different  dMRI  approaches. 

Diffusion  modeling  frameworks  may  be  roughly  categorized  into 
two  main  categories:  (1)  physical  or  “signal-driven”  diffusion  models, 
which  aim  at  measuring  and  characterizing  of  probabilistic  water  dis¬ 
placement  profile  (a.k.a.  the  diffusion  propagator);  and  (2)  biological  or 
“microstructure-driven”  models,  which  aim  to  provide  a  more  direct 
assessment  of  tissue  compartments  and  their  biological  attributes  by 
directly  incorporating  a  priori  biological  information  (assumptions)  into 
the  model.  Each  of  these  approaches  provides  a  set  of  scalar  outcome 
metrics  that  are  based  on  the  fitted  parameters  of  the  model,  and  the 
maps  of  these  scalars  are  used  to  visualize  and  quantify  differences 
across  brain  regions  or  to  identify  abnormalities  in  experimental 
models. 

The  first  class  of  techniques— physical  diffusion  modeling— is  the 
most  prevalent  and  includes  DTI  (Basser,  Mattiello,  &  LeBihan,  1994), 
diffusion  kurtosis  imaging  (DKI)  (Jensen,  Helpern,  Ramani,  Lu,  & 
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Kaczynski,  2005;  Tabesh,  Jensen,  Ardekani,  &  Helpern,  2011),  and  diffu¬ 
sion  propagator  representation  models  (Callaghan,  Eccles,  &  Xia,  1988; 
Callaghan,  MacGowan,  &  Packer,  1990).  These  models  directly  relate 
the  MRI  signal  for  each  DWI  to  the  diffusion  weighting  (i.e.,  b-value  in 
seconds  per  millimeter  squared)  and  the  orientation  of  the  applied  diffu¬ 
sion  sensitizing  gradients,  which  can  be  related  to  the  q-vector  (1/mm) 
in  reciprocal  displacement  space.  The  simplest  dMRI  approach  uses  as 
few  as  two  DWI  measurements  to  fit  a  single  exponential  decay  to  esti¬ 
mate  the  apparent  diffusion  coefficient  (mm2/s)  (Eccles,  Callaghan,  & 
Jenner,  1988)  of  water  present  in  the  tissue  in  the  direction  of  the 
applied  diffusion  weighting.  DTI,  the  most  commonly  used  diffusion 
model,  is  an  extension  of  this  approach  to  model  water  diffusion  in 
three  dimensions  by  the  diffusion  tensor  (Basser  et  al.,  1994).  The 
eigenvalues  (A.1;  X2,  and  X3)  of  the  diffusion  tensor  can  be  used  to  repre¬ 
sent  the  axes  of  a  diffusion  ellipsoid  related  to  the  preferred  orientation 
and  magnitude  of  water  diffusion  and  to  provide  a  set  of  scalar  metrics 
based  on  mathematical  arrangements  of  the  eigenvalues  that  report  a 
particular  feature  of  the  diffusion  tensor.  The  trace  of  the  diffusion  ten¬ 
sor  (TR  =  X1  +  X2  +  ^3)  or  mean  diffusivity  (MD  =  TR/3)  and  the  frac¬ 
tional  anisotropy  (FA)  (Pierpaoli  &  Basser,  1996)  are  the  most 
commonly  reported  DTI  scalar  metrics,  but  others,  including  axial 
(Dax  =  />Li)  and  radial  (Drad  =  [X2+XJi/2)  diffusivity  and  Westin  linear, 
planar,  and  spherical  anisotropy  (WL,  WP,  and  WS)  (Westin  et  al., 
2002),  can  also  be  applied  to  provide  more  specific  geometric  informa¬ 
tion  about  water  diffusion  derived  from  the  diffusion  tensor. 

While  DTI  already  offers  a  quantitative  and  sensitive  way  to  detect 
and  evaluate  changes  in  water  diffusion  that  arise  from  alterations  in 
the  tissue  environment,  the  diffusion  tensor  approximation  is  a  simpli¬ 
fied  representation  that  cannot  convey  higher-order  signal  behavior 
(i.e.,  non-Gaussianity)  or  complex  orientational  information  (i.e.,  multi¬ 
fiber  geometry).  To  address  these  limitations,  numerous  models  have 
been  developed  with  the  goal  of  quantifying  the  diffusion  propagator,  P 
(r) — the  3D  probability  distribution  for  water  displacement  (Callaghan 
et  al.,  1990).  The  earliest  propagator  measurement  approaches  relied 
directly  on  the  Fourier  relationship  between  P(r)  and  the  MRI  signal  of 
DWI  data  collected  in  q-space  E(q)  (Callaghan  et  al.,  1988).  Q-space 
imaging  methods  of  this  type  are  model-free  and  simply  rely  on  taking 
a  mathematical  transform  of  E(q)  at  each  voxel.  These  include  diffusion 
spectrum  imaging  and  Q-ball  imaging  (Tuch,  2004;  Tuch,  Reese,  Wie- 
gell,  &  Wedeen,  2003),  which  produce  orientation  distribution  func¬ 
tions  (ODFs)  that  allow  visualization  of  different  fiber  orientations  in 
the  same  voxel. 

More  recently,  model-dependent  approaches  to  propagator  repre¬ 
sentation  have  been  proposed  that  report  information  about  water  dif¬ 
fusion  that  is  “non-Gaussian."  For  example,  DKI  (Glenn,  Helpern, 
Tabesh,  &  Jensen,  2015;  Jensen  et  al.,  2005;  Tabesh  et  al.,  2011) 
expands  the  DTI  model  with  an  additional  term  (Table  2)  that  includes 
the  kurtosis  tensor,  and  several  metrics  may  be  calculated  from  this  new 
tensor  including  the  mean,  radial,  and  axial  kurtosis  (MK,  RK,  and  AK) 
as  well  as  the  kurtosis  FA  (KFA).  In  addition  to  new  DKI  metrics  that 
probe  non-Gaussian  information  in  the  dMRI  signal,  the  conventional 
DTI  metrics  that  are  calculated  using  this  model  are  more  stable  than 
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for  DTI  alone  (Hutchinson  et  al.,  2017;  Veraart  et  al..  2011),  and  the 
DWI  sampling  for  DKI  is  modest  compared  with  higher-order  models 
and  therefore  requires  relatively  shorter  scan  times.  However,  DKI  is 
highly  dependent  on  the  DWI  sampling  scheme,  potentially  leading  to 
inconsistent  DKI  values  across  studies.  Several  studies  of  CCI  in  the  rat 
have  suggested  that  DKI  markers,  especially  increases  in  MK,  may  be 
more  sensitive  to  particular  abnormalities  such  as  the  lesion  boundary 
(Jiang  et  al.,  2011)  and  the  temporal  evolution  of  abnormalities,  which 
is  distinct  from  the  pattern  of  DTI  changes  (Zhuo  et  al.,  2012).  How¬ 
ever,  DKI  has  also  been  observed  to  be  less  sensitive  to  abnormalities 
in  a  repetitive  mild  TBI  model  (Yu  et  al.,  2017),  suggesting  that  in  some 
cases  the  benefit  of  increased  specificity  is  offset  by  reduced 
sensitivity. 

While  DKI  assumes  a  particular  signal  behavior  for  all  non-Gaussian 
diffusion,  other  higher-order  models  have  been  developed  that  are 
more  flexible  and  extensive  for  describing  non-Gaussian  diffusion.  For 
example,  mean  apparent  propagator  (MAP)-MRI  (Avram  et  al.,  2016; 
Ozarslan  et  al.,  2013)  and  ensemble  average  propagator  approaches 
(Cheng,  Jiang,  &  Deriche,  2012;  Hosseinbor,  Chung,  Wu,  &  Alexander, 
2013)  use  carefully  selected  functional  series  expansion  to  fit  DWI  data 
sampled  over  a  wide  range  of  b-values  with  many  orientations.  The 
coefficients  of  these  models  provide  an  approximation  of  the  diffusion 
propagator  that  may  be  visualized  as  ODFs  or  generate  scalar  maps  of 
relevant  propagator  metrics.  For  MAP-MRI,  these  include  zero- 
displacement  probabilities  (return  to  origin,  axis,  and  plane  probability, 
or  RTOP,  RTAP,  and  RTPP,  respectively),  non-Gaussianity  (NG),  and 
propagator  anisotropy  (PA).  While  the  greater  complexity  of  techniques 
that  measure  the  propagator  provides  more  detailed  information  about 
water  diffusion  that  may  lead  to  improved  specificity,  these  models 
require  considerably  more  DWI  data  with  greater  diffusion  weighting. 

While  physical  diffusion  models  present  scalar  metrics  that  are 
related  to  the  extent  and  preferred  orientation  of  water  diffusion,  bio¬ 
physical  models  incorporate  a  priori  assumptions  about  the  tissue  envi¬ 
ronment  to  estimate  features  such  as  cellular  fraction  volumes,  axon 
and  tract  dimensions,  or  neurite  density  and  coherence.  These 
approaches  are  appealing  for  identification  of  abnormalities  in  biologi¬ 
cal  terms;  however,  caution  is  warranted  not  only  to  consider  the  valid¬ 
ity  of  model  assumptions  across  different  healthy  tissue  types  but  also 
to  consider  how  the  biophysical  model  does  or  does  not  accommodate 
pathological  processes.  For  example,  tractography  approaches  (Tour- 
nier,  Calamante,  &  Connelly,  2012;  Tournier,  Calamante,  Gadian,  & 
Connelly,  2004),  which  represent  white  matter  “tracts”  using  DTI  or 
ODF  information,  may  misrepresent  pathological  changes  that  increase 
tissue  organization,  such  as  organized  gliosis  as  a  false  increase  in  tract 
density  or  length,  or  other  pathological  changes  such  as  gliosis  may 
decrease  white  matter  anisotropy  even  if  the  tracts  are  still  present, 
which  may  lead  to  a  misleading  loss  of  tract  representation  (see  Figure 
3  and  references;  Budde  et  al.,  2011).  This  said,  tractography  can  be 
used  effectively  to  illustrate  white  matter  differences  or  to  delineate 
known  tract  anatomy  when  interpretation  is  made  carefully  with  refer¬ 
ence  to  known  neuroanatomy  and/or  accompanied  by  histological  or 
functional  validation. 
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Other  biophysical  methods  have  been  developed  to  provide  scalar 
metrics  that  report  local  information  about  the  tissue  environment.  For 
example,  the  composite  hindered  and  restricted  model  of  diffusion 
(CHARMED)  (Assaf  &  Basser,  2005)  and,  later,  AxCaliber  (Assaf, 
Blumenfeld-Katzir,  Yovel,  &  Basser,  2008)  were  developed  to  estimate 
white  matter  fiber  orientation  and  axon  diameter  in  part  to  improve 
the  sensitivity  of  dMRI  to  detect  abnormalities  in  disorders  affecting 
the  white  matter  and  may  have  particular  relevance  for  evaluating  axo¬ 
nal  damage  and  myelination  changes.  Cellular  compartment  models 
based  on  multiexponential  dMRI  modeling  or  more  advanced 
approaches  have  also  been  devised  to  estimate  the  intra-  and  extracel¬ 
lular  contribution  to  the  dMRI  signal  and  to  further  model  axonal 
dimensions  (Jelescu,  Veraart,  Fieremans,  &  Novikov,  2016;  Jelescu, 
Zurek,  Winters,  &  Veraart,  2016;  Kunz  et  al.,  2014;  Panagiotaki  et  al., 
2012).  More  recently,  several  models  have  extended  compartmental 
modeling  to  include  features  of  “white  matter  tract  integrity”  (WMTI; 
Fieremans,  Jensen,  &  Helpern,  2011)  or  “neurite  orientation  dispersion 
and  density”  (NODDI;  Tariq,  Schneider,  Alexander,  Gandini  Wheeler- 
Kingshott,  &  Zhang,  2016;  Zhang,  Schneider,  Wheeler-Kingshott,  & 
Alexander,  2012),  which  may  be  advantageous  for  the  detection  of 
pathological  or  regenerative  changes  to  specific  neuroanatomical  struc¬ 
tures  (i.e.,  white  matter  or  neurites).  While  initial  work  using  biophysical 
modeling  has  demonstrated  correlations  between  histological  and 
dMRI-derived  values  for  white  matter  axonal  fraction  and  diameter  fol¬ 
lowing  CCI  in  rats  (Wang  et  al.,  2013),  much  remains  to  be  understood 
about  the  ability  of  biophysical  models  to  accurately  report  tissue 
parameters.  Furthermore,  biophysical  models  have  been  shown  to 
depend  on  specific  fixed  model  parameters  (e.g.,  compartmental  diffu- 
sivity  values),  DWI  sampling,  and  image  quality  (Jelescu,  Veraart,  et  al., 
2016),  such  that  extreme  care  should  be  given  to  application  of  these 
models. 

4  I  A  PRACTICAL  AND  CONCEPTUAL 
PERSPECTIVE  FOR  THE  USE  OF  ADVANCED 
DMRI  METHODS  IN  TBI  RESEARCH 

The  past  several  decades  have  seen  considerable  growth  in  our  ability 
to  characterize  TBI-related  brain  changes  by  incorporating  dMRI  into 
studies  of  experimental  brain  injury  models.  From  early  work  to  define 
diffusivity  changes  related  to  pathophysiology  following  TBI  to  more 
recent  endeavors  to  relate  DTI  abnormalities  with  a  range  of  cellular 
alterations,  dMRI  has  emerged  as  a  promising  tool  for  providing  micro¬ 
scale  information  about  brain  abnormalities  that  is  quantitative,  whole- 
brain,  and  noninvasive. 

The  future  of  these  efforts  will  be  better  enabled  if  dMRI 
approaches  can  be  used  in  combination  with  one  another  and  along¬ 
side  other  MRI  modalities  as  a  set  of  specialized  tools  to  detect  specific 
pathology,  predict  outcomes,  or  target  treatments.  However,  for  this  to 
become  possible  we  must  first  systematically  understand  the  similar¬ 
ities  and  differences  across  dMRI  models  including  their  relative  advan¬ 
tages  and  limitations.  A  recent  study  comparing  the  dependence  of 
scalar  metrics  from  DTI,  DKI,  MAP-MRI,  and  NODDI  in  normal  brain 


tissue  identified  several  practical  implications  of  moving  to  models  with 
greater  complexity  or  biological  assumptions  (Hutchinson  et  al.,  2017). 
In  particular,  DKI  metrics  were  found  to  be  the  most  vulnerable  to 
DWI  sampling  and  image  quality,  while  MAP-MRI  required  the  greatest 
DWI  sampling  and  was  dependent  on  the  initial  fitting  of  the  diffusion 
tensor.  The  NODDI  model  was  found  to  have  stability  over  a  range  of 
DWI  sampling  and  image  quality,  but  to  depend  greatly  on  the  selec¬ 
tion  of  fixed  model  parameters  (e.g.,  interstitial  diffusivity).  Greater 
understanding  of  the  influence  of  experimental  parameters  on  dMRI 
outcome  measures  as  well  as  identifying  redundancy  and  novelty  of 
information  across  models  is  an  active  area  of  imaging  research  and 
will  aid  in  the  effective  use  of  these  methods. 

Despite  the  challenges  that  accompany  more  complex  diffusion 
modeling,  there  is  also  promise  in  their  ability  to  better  characterize 
posttraumatic  brain  abnormalities  than  existing  methods.  This  was 
recently  demonstrated  in  a  study  combining  DTI  and  DKI  metrics  to 
evaluate  markers  of  pathology  after  experimental  TBI  in  rats  (Zhuo 
et  al.,  2012).  While  cortical  MD  was  decreased  2  hr  after  injury  and 
then  increased  1  week  later,  the  MK  was  consistently  elevated  at  both 
times  and  furthermore  corresponded  to  astrocytosis.  While  more  work 
remains  to  define  whether  and  how  DKI  and  other  advanced  methods 
provide  additional  or  more  specific  information  that  corresponds  to 
individual  or  organized  pathogenesis,  this  study  and  others  employing 
advanced  dMRI  methods  to  evaluate  experimental  TBI  (Davoodi-Bojd 
et  al.,  2014;  Fozouni  et  al.,  2013;  Wang  et  al.,  2013;  Yu  et  al.,  2017) 
have  begun  to  build  an  important  foundation  to  understand  and  apply 
advanced  diffusion  imaging  in  the  context  of  TBI  research. 

The  combination  of  a  modern  dMRI  toolkit  with  other  MRI  modal¬ 
ities  and  advances  in  the  understanding  of  neurobiological  responses 
to  TBI  has  the  potential  to  improve  the  spatiotemporal  characterization 
of  TBI  across  various  experimental  models  and  provide  outcome  meas¬ 
ures  that  may  be  directly  translated  as  imaging  markers  in  human  stud¬ 
ies  or  used  in  experimental  models  for  the  development  of  therapies. 
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